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We present a Hartree-Fock study that incorporates the ef- 
fects of Landau level mixing and screening due to filled levels 
into the computation of energies and states of quasiparticles 
in quantum Hall ferromagnets. We use it to construct a phase 
diagram for skyrmion stability as a function of magnetic field 
and Zeeman coupling strengths. We find that Landau level 
mixing tends to favor spin-polarized quasiparticles, while fi- 
nite thickness corrections favor skyrmions. Our studies show 
that skyrmion stability in high Landau levels is very sensitive 
to the way in which electron-electron interactions are modified 
by finite thickness, and indicate that it is crucial to use mod- 
els with realistic short distance behavior to get qualitatively 
correct results. We find that recent experimental evidence for 
skyrmions in higher Landau levels cannot be explained within 
our model. 
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I. INTRODUCTION 

It has been recognized that exotic magnetic excitations 
known as skyrmions may existuu in a two-dimensional 
electron gas in a strong homogeneous magnetic field 
(quantum Hall system) near spin polarized groundstates. 
These are excitations of a two-dimensional spontaneous 
ferromagnet, the physics of which is relevant to this sys- 
tem (despite the presence of a strong magnetic field), 
because of the small Lande g-factor in GaAs systems 
(where most experiments take place), which makes the 
Zeeman coupling very small compared to other energy 
scales (Coulomb interaction, cyclotron energy) in the 
problem. Skyrmions are spin configurations with a non- 
trivial winding number (Pontryagin index). They were 
first discussed in the context of four-dimensional field 
theories,QEl and were later recognized as states occur- 
ring in the non-linear sigma model description of two- 
dimensional ferromagnets. □ For filling factors v = N/N<j, 
close to one (N is the number of electrons and N^, the 
number of magnetic flux quanta penetrating the system), 
these turn out to be the lowest energy quasiparticles un- 
der typical experimental circumstances. Skyrmions can 
thus be introduced into the groundstate by adding or 
removing charge from the system.0 

Experimentally, the case for the existence of skyrmions 
in a system close to v — 1 is quite strong. NMR exper- 
iments show a degrading of the spin polarization with 
deviation of filling factor from ondj that is in .teinark- 
ably good agreement with Hartree-Fock theory BlI The. 
quasiparticle spin measured in transport experimentsBI 



are also reasonably well accounted for by Hartree-Fock 
calculations. a Electromagnetic absorption experimentsE2l 
further support that doping away from v = 1 injects 
skyrmions into the system. 

In weaker magnetic fields, near filling factor v = 3, 
early experimentsO suggested that spin-polarized quasi- 
particles are lower in energy than skyrmions, so the ef- 
fects seen near v = 1 would not be present at higher 
filling. This is consistent with calculations of skyrmion 
energiesjiear v = 3 that include finite thickness correc- 
tions,!!^ c3 which indicate that skyrmions will be present 
only at much smaller Zeeman couplings than realized in 
typical experiments. The size of the skyrmion may be 
quantified by a number K, the difference in the spin com- 
ponent S z between the skyrmion and the spin-polarized 
quasiparticle. Because of the necessarily small Zeeman 
coupling, stable skyrmions close to v — 3 have large val- 
ues of K. They also become unstable with respect to 
spin-polarized quasiparticles at a finite value of K. (For 
v = 1, K — *■ as the Zeeman coupling reaches the maxi- 
mum value for which the system supports skyrmions; i.e., 
the skyrmion state smoothly goes into the spin-polarized 
quasiparticle state.) For a two-dimensional electron gas 
(2DEG) with width of about 2£, where I = ^JheJeB is 
the magnetic length, the minimum K expectedlj is ap- 
proximately 4. J—. 

Recently however, NMR experimentsEJ have uncovered 
evidence that some anomalous degrading of spin polar- 
ization does occur as one dopes away from v = 3 at rela- 
tively high Zeeman couplings. These experiments further 
indicate that the number of overturned spins per quasi- 
particle is quite—small, K ~ 1. The simple models usu- 
ally consideredtlrEj are inconsistent with this, and one 
is naturally led to inquire as to what other ingredients 
might change the critical Zeeman coupling and small- 
est K observable near v = 3. Two possible answers are 
Landau level mixing and screening by filled Landau lev- 
els. It should be noted that these are not distinct effects: 
screening by filled Landau levels occurs because they may 
admix high (unoccupied) Landau levels to smooth fluctu- 
ations due to external potentials and/or inhomogeneous 
electron densities in partially filled levels. Conversely, 
the states which may be used for Landau level mixing in 
a partially filled level are limited to those that are not 
occupied by electrons in other levels, due to Pauli ex- 
clusion. Thus, a correct treatment of either screening or 
Landau level mixing near v — 3 must include both these 
effects. 

In this work, we present a method by which these 
may be incorporated into the Hartree-Fock description 



of skyrmion states. 

Our principal conclusions may be summarized as fol- 
lows: (i) For skyrmions near v = 1, Landau level mix- 
ing tends to lower the quasiparticle energy, although not 
enough to quantitatively explain the activation energies 
seen in experiment!! Introduction of a finite width of one 
magnetic length lowers the energy of the skyrmion by 
approximately 40%, and inclusion of Landau level mix- 
ing lowers the energy by approximately another 20% for 
Tiuj c /(e 2 / k£) « 0.6, where uj c = eB/m*c is the cyclotron 
frequency of the electrons, and k is the dielectric con- 
stant of the host crystal. The resulting quasiparticle 
gap is approximately|-.a factor of 2 larger than what is 
found in experiment Ij This result agrees qualitatively 
with that of another study of Landau level mixing ef- 
fects on v = 1 skyrmions 113 (ii) For v = 1, Landau level 
mixing tends to suppress quasihole-like skyrmions (i.e., 
lowering the maximum Zeeman coupling for which they 
are stable), while enhancing the stability of quasielec- 
tron (anti-)skyrmions. (iii) For v — 3 and higher, we find 
that a sufficiently realistic model of the effective electron- 
electron interaction, as modified by the finite thickness of 
the electronic wavefunctions, is necessary to obtain reli- 
able results—Use of a simple potential due to Zhang and 
Das SarmaiLa grossly overestimates rUae-stability of v = 3 
skyrmions; more realistic potentialail'tHl allow skyrmions 
only for very small Zeeman couplings, (iv) Screening and 
Landau level mixing for v — 3 and v = 5 tend to lower the 
energy of spin-polarized quasiparticles more than that of 
skyrmions, making the latter even less stable, (v) The 
results of Ref. [lj cannot be understood solely on the 
basis of Hartree-Fock states for skyrmions. 

The remainder of this article is organized as follows. 
In Section II below, we discuss the method used to allow 
screening and Landau level mixing to be included in our 
calculations. Section III gives details of our results, and 
we conclude with a summary of our findings in Section 
IV. 



II. HARTREE-FOCK METHOD WITH LANDAU 
LEVEL MIXING 

Most previous Hartree-Fock studies of skyrmions have 
relied on Landau level representation of the single par- 
ticle states .B We choose instead to construct the wave- 
functions in real space. This enables us to include in the 
model Landau level mixing occurring in weak magnetic 
fields, without having to expand over the large number 
of Landau levels necessary in the former approach. We 
thus trade the calculational convenience of working with 
the functions given in closed analytic form (Landau lev- 
els) for a closer description of the single-particle states 
by representing them on a real-space grid. 

In this calculation we aim to model all the participat- 
ing particles, including the ones in the filled levels. Our 
Hartree-Fock wavefunction is a Slater determinant com- 
posed of single-particle states which have L z ± S z as a 



good quantum numberi but whose radial form is to be 
determined self-consistently: 



^skyrmion) = Y\l\ m \^) 

fim(r) 
g lm (r)e ±M 



(1) 



Here r and 9 are polar coordinates, m is the angular 
momentum quantum number, and i labels different states 
of the same m. The sign ± corresponds to two families 
of solutions, + for antiskyrmion (or quasielectron spin 
structured solution) and — for skyrmion (quasihole). 

In very strong magnetic fields, the functions f(r) and 
g(r) take the form expected for Landau level states. 
When the strength of the magnetic field is lowered to 
bring the ratio of cyclotron and Coulomb energy scales 
close to 1 , the form of the radial part relaxes toward 
some modified form, as dictated by the interactions in 
the system. 

Using the trial form of the wavefunction (jl]) , the many- 
body Schrodinger equation with the Hamiltonian 

ff=-L / d 2 rV*t (f) |^V--!| 2 * CT (r) 
2m J ' % c 

(7 



2^ 



d z rd 2 r'^l(r)m a ,{r)v{r- r')^,(r')* CT (r) (2) 



(where a denotes spin and v(f — r 1 ) the Coulomb interac- 
tion), upon variation with respect to the functions / and 
g, gives a system of mean-field single-particle equations: 

D m (r)fim(r) - ^g^BBf lm (r) 

+ / r'dr'V H (r,r')[p(r')~p ]f lm (r) 
Jo 



r'dr'^V£L m ,(ry)p%(r',r)f im (r') 



r> dr' ]T C-w (r, r>£ (r', r)g im±1 {r') 



(3) 



together with the analogous equation for the function 
g(r). Here is a dictionary of the notation accompanying 
Eq. |^: the operator D rn is given by 



h 1 d d 

Dm = — T r ~r 

2m* r dr dr 



flLJr 



(m*) 2 Lj 2 c r' 
4 2m* 



(4) 



With LO r 



B the magnitude of the external magnetic 
field, and m* the effective mass of the electron, g is the 



Lande g-factor, and /is is the Bohr magneton, p's denote 
generalized densities: 

P^(r',r) = f^(r')fm'(r) 
P™>(r',r) = g* m , ±1 {r')g m , ±l {r) 
P i Mr',r)=g* m , ±1 {r')f„Ar) 

p(r) = J2^M r ' r )+p l M r ^ ( 5 ) 



and po is the uniform background density. V H and V ex 
are the following integrals of the Coulomb potential over 
the azimuthal variable: 

r2lT p2tt 

V H (r,r')= del dO'v{r-P) 
Jo Jo 

ex / , m I d6'e l ( m - m 'W~ 9 ">v(r-P), (6) 

^0 



V%L m ,(r,r')= d9 



and, finally, ti stands for the single-particle Hartree-Fock 
energy. 

The finite width of the sample is modeled jawing the 
form of the in-plane potential due to Cooper :til we re- 
place the Coulomb interaction v(f ' — r') in the Eq. @ 
by 



vc(r) 



dzdz' - 



-(z' 2 +z 2 )/2w 2 



1 



2ttw 



^r 2 + {z- z') 



A2 



(7) 



The symbol w denotes the width of system in the direc- 
tion perpendicular do the plane of the system, and z, z' 
are coordinates in that direction. 

To handle the boundaries of the system, we assume the 
electron states with angular momentum m > m max have 
the ferromagnetic groundstate form (i.e. Landau levels 
with well defined spin). The states with m < m max are 
explicitly included in the calculation. For K not too large 
we find it is sufficient to allow variations of the states 
with m of up to 30 for v = 3 and up to 50 for v = 5. 
In practice, including boundary electrons from the states 
with m between 31 and 100 (y = 3) and between 51 and 
120 (v = 5) describes the effect of the system edge with 
precision matching the rest of the calculation. 

Understanding Eq. (||) as a system of coupled eigen- 
problems, we look for the self-consistent single particle 
solutions. (Discretization will turn each eigenequation 
into a matrix diagonalization problem which can be han- 
dled using standard methods.) The results we thus ob- 
tain will be largely presented as comparisons between 
energies of the spin-polarized quasiparticle and energies 
of the corresponding skyrmion. 

To assess the energy of the skyrmion in the region of 
parameters where it is not stable, we add to the Hamil- 
tonian a term of the form X(S Z — So) 2 , S z being the spin 
operator, and A a tunable parameter. This term favors a 
state with total spin So, but is insensitive to the detailed 
form of the wavefunctions. This allows the variational 
scheme to pick out the lowest energy Slater determinant 



of the form given in Eq. (jl|) within the space of states with 
the same fixed value of K. 



III. RESULTS 

Based on the calculation described in the previous sec- 
tion, we present some of the results the method allows 
us to obtain; we focus mainly on the singly charged ex- 
citations in the first three Landau levels. Consideration 
of higher Landau levels is also possible, but computation 
of the potential lookup tables becomes prohibitively ex- 
pensive, and, as the results so far indicate, leads to no 
new insight. In the following we shall take the unit of 
energy to be e 2 / k£, and the unitless Zeeman splitting to 
be g = 9 JiJ^ , where e is the electron charge, k is the 
dielectric constant of the host material, I is the magnetic 
length in the field B and /is is the Bohr magneton. 



A. Skyrmion vs. polarized quasiparticle 

In Fig. we show the energy difference between the 
spin-polarized quasiparticle and the skyrmion of size K, 
Vk — Vk=o- This quantity is a pure interaction energy 
(i.e. Zeeman energy is not included), and represents the 
energy gained or lost in deforming a spin-polarized quasi- 
particle into a skyrmion when Zeeman coupling is absent. 
Of particular importance is the slope (negative slopes in- 
dicate that skyrmion is stable for some value of g) and 
the curvature (concave curves will support small-sized 
skyrmions). For concave curves the largest Zeeman split- 
ting that supports .skyrmions is the negative of the initial 
slope of the curve.E 2 ] 

For large cyclotron energies our results are essentially 
identicakto those obtained using the single Landau level 
method.™ Note the quasielectron and quasihole exci- 
tations are precisely degenerate in this case, due to 
particle- hole symmetry. For smaller values of %u) c , the 
two curves split; surprisingly, the quasihole skyrmion is 
suppressed by Landau level mixing, whereas the quasi- 
electron skyrmion is enhanced. (The former result is in 
agreement with Ref. Oft. 
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FIG. 1. Energy difference between skyrmion and polar- 
ized quasihole at v = 1, for well width of 1.73£ and with 
g — * 0. The full lines correspond to the quasihole and the 
dashed ones to the quasielectron case. For high cyclotron en- 
ergies (100. Oe 2 / k£ ; the lower full curve) the results for the 
two types of excitations overlap. For low cyclotron energies 
upper full line and non-overlapping dashed line) 
there is significant difference in the behavior of the two, which 
reflects itself in the phase diagram, Fig. ^[ 

The energy gaps (Fig. ||) which result from creation of 
skyrmion-antiskyrmion pairs when Landau level mixing 
and finite thickness corrections are included are consid- 
erably smaller than what is found for two-dimensional 
layers and no mixing.0 However, the resulting energies 
are still almost-a factor of two larger than what is found 
in experiment .□ The discrepancy is likely to be due to 
disorder. 




FIG. 2. Skyrmion-antiskyrmion pair excitation gap for 
two different well widths : 1.73^ (full line) and 1.01 (dashed 
line). Upper lines correspond to the cyclotron energy of 
100.0e 2 /fi;f and lower ones to 0.8e 2 /«;f . 



Figs. (||) and (||) present analogous results for v = 3 
and v — 5. Note the considerably smaller energy scales in 
these figures, indicating that skvMiiians can only be sta- 
ble (if ever) for small values of gO'EHl It is apparent that 
the introduction of Landau level mixing and screening 
destabilizes the skyrmion. Evidently, spin-polarized par- 
ticles are better able to take advantage of the admixture 
of higher Landau levels than skyrmions. 
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FIG. 3. Energy difference between skyrmion and polarized 
quasihole at v = 3; the well width is 1.73^ ; g — * 0. The 
full lines correspond to the quasihole and the dashed ones to 
the quasielectron case. The cyclotron energy range is between 
lOO.Oe 2 / k£ for the lowest line and 0.8e 2 j kJL for the uppermost, 
each cyclotron energy being 50% smaller than the previous 
higher one. 
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FIG. 4. Energy difference between skyrmion and polarized 
quasihole at v — 5; the well width is 1.731 ; g — > 0. The full 
lines correspond to the quasihole and the dashed ones to the 
quasielectron case. The cyclotron energies are the same as in 
Fig. | 



B. Effect of Finite Well Width 

Quasiparticle energies depend on the waUj-syidth, as 
illustrated in Figs. || and [| As expect edJIirEa we find 
that for narrower wells the difference in energy is less 
favorable for the skyrmion (Fig. ^j). Note that the width 
used in. Fig. |^ is close to an experimentally reported 
value o 



vzds(r- f") 



■>J\r — F| 2 



(8) 
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FIG. 5. Energy difference between skyrmion and polarized 
quasihole at v — 3 and the well width of 1.73^ ; g — > 0. The 
cyclotron energies are the same as in Fig. H. 
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K (minority spin) 

FIG. 6. Energy difference between skyrmion and polarized 
quasihole at v — 3 and the well width of 1.0£ ; g — > 0. The 
cyclotron energies are the same as in Fig. 0. 



which is commonly used in studying quantum Hall sys- 
tems (see for example Refs. [l3]and[l(]) gives substantially 
different results then those presented above (Fig. H). 
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FIG. 7. Difference in energy between the skyrmionic and 
polarized quasihole solutions using Zhang-Das Sarma poten- 
tial. The well width is 1.73£ and v — 3. The cyclotron energies 
are the same as in Fig. |3| 

The principal difference between vc(r) and vzds{ r ) is 
the behavior at small r; the former diverges logarithmi- 
cally, whereas the latter is regular. Other divergent po- 
tentials give results consistent with Figs. || and [|; it is 
likely that the oversimplified behavior at short distances 
is responsible for the poor performance of vzdS ( r ) m this 
problem. 



C. Phase Diagram 

Based on results in Figs. [|, S, and |J we can con- 
struct the phase diagrams of skyrmion stability for the 
filling factors of v = 1, 3, and 5. Large uj c and small g c is 
the region favoring the spin structured excitations. We 
see that the region "shrinks" as one moves to the higher 
filling factors. Also, according to this calculation the 
breaking of symmetry between the quasihole and quasi- 
electron excitations upon lowering the cyclotron energy 
is quite spectacular in the lowest Landau level, whereas 
it plays no significant role in the higher ones. 



It is worth remarking at this juncture that a reasonably 
realistic model of the electron-electron interaction with 
finite sample thickness corrections is needed to obtain 
qualitatively correct results. Fig. fj] shows that thejise 
of a simpler model potential (Zhang and Das SarmaEj) 




FIG. 8. Phase diagram of skyrmion stability. The dashed 
lines correspond to the the spin structured quasielectron case, 
whereas the full lines correspond the quasihole one. Top lines 
correspond to the case of v = 1, middle to v = 3 and bottom 
to v — 5. The scale on the abscissa is logarithmic. The width 
is w = 1.731 

In the paper by Song et al. the reported excitation 
at the parameter values of w = l.73£, uj c = 0.56, and 
g = 1.72 x 10~ 2 falls well outside the boundary expected 
from the Hartree-Fock calculation. 



D. Effect of Impurities 

It is tempting to speculate that inclusion of impurity 
effects can stabilize the skyrmions at values of g bigger 
than allowed in a pure sample. To test this idea we can 
include a simple model of an impurity in our calculation: 
a point charge (impurity) at a distance zq above the cen- 
tral plane of the system. It is replaced by an effective 
non-uniform charge density in the plane producing the 
same potential, 

The effective density can be found to be 

1 Zq 

^ (r)= Mz 2 + i?2)3/ 2 ( 10 ) 

for zq > 0. Results with and without such an impurity 
are illustrated in Fig. ^. As may be seen, the impu- 
rity favors spin-polarized quasiparticle over skyrmion. A 
similar result is expected for a short-range (e.g. delta- 
function) impurity potential. Apparently the simplest 
models of disorder are not likely to explain the results 
of Ref. 0~3. It is probable that more complicated im- 
purities (e7g. multiply charged or magnetic ones) could 
stabilize the small-spin skyrmions at v — 3. However, in 
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FIG. 9. Comparison of skyrmion stability in a system 
without an impurity (dashed line) and the one with impurity 
at the distance of 2 magnetic lengths from the sample for 
v — 1 case. Cyclotron energy is 100. Oe 2 //^ in both cases. 
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1.1 2.2 3.3 4.4 

K (minority spin) 

FIG. 10. The same as Fig. | for v = 3 case. 

IV. CONCLUSION 

In this paper we have presented a real-space method 
for computing Hartree-Fock states and energies of two- 
dimensional systems in magnetic fields, appropriate for 
systems with circular symmetry in which Landau level 
mixing may be important. The method was applied to 
compute the effects of Landau level mixing and screening 
on skyrmion states. It was found that in most cases these 
tend to destabilize skyrmions, with a notable exception 



occurring for the case of the quasielectron (antiskyrmion) 
around v = 1. The calculations indicate that Hartree- 
Fock states cannot account for the results of Ref. (re- 
porting skyrmions at v = 3). This is in agreement with 
earlier studies where Landau level mixing and screening 
were not included. 
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